Linking numbers for self-avoiding loops and percolation: 
application to the spin quantum Hall transition 
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Non-local twist operators are introduced for the O(n) and Q-state Potts models in two di- 
mensions which, in the limits n — > (resp. Q — > 1) count the numbers of self-avoiding loops 
(resp. percolation clusters) surrounding a given point. Their scaling dimensions are computed ex- 
actly. This yields many results, for example the distribution of the number of percolation clusters 
which must be crossed to connect a given point to an infinitely distant boundary. Its mean behaves 
as (l/3v37r)| m(p c — p)\ as p — > p c — . These twist operators correspond to (r, s) = (1,2) in the 
Kac classification of conformal field theory, so that their higher-order correlation functions, which 
describe linking numbers around multiple points, may be computed exactly. As an application we 
compute the exact value yS/2 for the dimensionless conductivity at the spin Hall transition, as well 
as the shape dependence of the mean conductance in an arbitrary simply connected geometry with 
two extended edge contacts. 

PACS numbers: 05.50.+q, 64.60.Ak, 73.40.Hm 

The conformal field theory/Coulomb gas approach to 
two-dimensional percolation and self-avoiding walk prob- 
lems has been extraordinarily fruitful 0-0] • In addition 
to values for many of the critical exponents, other univer- 
sal scaling functions such as percolation crossing proba- 
bilities have been obtained exactly J|||. In this Letter 
a set of correlation functions of non-local operators is 
introduced, which describe topological properties of per- 
colation clusters and self- avoiding walks, and count the 
number of clusters or loops which must be crossed in or- 
der to connect two or more points. 

It turns out that these exponents may easily be com- 
puted using standard Coulomb gas methods, for general 
n or Q in the O(n) or Q-state Potts model respectively. 
In the limits n — > or Q — > 1, corresponding to self- 
avoiding walks or to percolation respectively, the scaling 
dimensions of these operators vanish, so that some of 
their correlations are trivial. However, it is their deriva- 
tives with respect to n or Q which give physical infor- 
mation, thereby giving rise to a variety of logarithmic 
behavior. The frequent occurrence of logarithmic cor- 
relations in such conformal field theories (CFTs) with 
vanishing central charge has recently been pointed out in 
several contexts 



In certain cases these topological operators may be rec- 
ognized as twist operators which have already been iden- 
tified in c = 1 theories || and the 3-state Potts model 
||. From the CFT point of view, they turn out to cor- 
respond to degenerate Virasoro representations labeled 
by (r, s) — (1,2) in the Kac classification. These bulk 
operators have not previously been identified for general 
n and Q. The fact that they are degenerate means that 
their higher-point correlations may be computed exactly. 

In addition to bulk operators of the above type it is also 
possible to define operators which count the number of 



loops or clusters surrounding a point near the boundary 
of a system. They may be used to compute the universal 
mean conductance at the spin quantum Hall transition, 
which has recently been shown to map onto the percola- 
tion problem . 

0(n) model. Let us recall the elements of the Coulomb 
gas approach n-component spins s(r) with s 2 (r) = 1 
are placed at the sites r of a lattice, with a nearest 
neighbor interaction. The partition function Tr Yl rr ' {^~^ 
ys(r) ■ s(r')J, when expanded in powers of y, gives a sum 
over self-avoiding loops with a factor y for each bond and 
n for each loop. Each loop may be replaced by a sum over 
its orientations, with a weight e x for each, with \ cho- 
sen so that the sum gives n = 2cos7rx. This gas of ori- 
ented loops is then mapped onto a height model with vari- 
ables 4>{R) G 7rZ on the dual lattice, such that on the dual 
bond RR' 4>(R) — 4>(R') = 0, ±1 according to whether the 
bond it crosses is empty or is occupied by a loop segment 
of one orientation or the other. At the critical fugac- 
ity y c this is supposed to renormalize onto a free field 
with reduced free energy functional (<?/(47r)) J(d<j)) 2 d 2 r, 
with respect to which the long-distance behavior of all 
correlations may be computed as long as phase factors 
associated with non-contractible loops are correctly ac- 
counted for. g may be fixed by a variety of methods: 
for the dilute regime of interest here, g = 1 + X- On & n 
infinitely long cylinder of perimeter L, in order to cor- 
rectly count loops which wrap around the cylinder with 
weight n it is necessary to insert factors e ±lx ^^ ±oc ' ) at ei- 
ther end: these modify the free energy per unit length to 
— (7r/6i)(l — (6/g)x 2 )i from which the value of the cen- 
tral charge c = 1 — 6(g — l) 2 /g follows. For self-avoiding 
walks, n = 0, x — hi an d g = §, so that c = 0. 

Now suppose that loops which wrap around the cylin- 
der are counted with a different weight n' = 2cos7rx'. 
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(A similar construction was used in Ref. (TT) to count the 
winding angle of open walks). The cylinder free energy 
per unit length will be modified by a term {ir/g){x! — X 2 )- 
In the plane, this may be interpreted as the insertion of 
a non-local operator a n i whose scaling dimension is 



x(n;n') = (l/2g(n))( X ' 2 -X 2 )- 



(1) 



so that the two-point function (av(ri)<j n '(r2)) decays 
like \r% — r2 \- 2x ( n > n ) a t criticality. The interpretation 
of this is as a defect line joining r\ and T2'. loops which 
cross this defect line an odd number of times, that is, 
which surround either but not both of the points Jl2[ , 
now carry a factor n' instead of n. A particular example 
is n' — —n, which is equivalent to making the identifi- 
cation s(r) — > — s(r) as the defect line is crossed. Such 
twist operators have been identified in other models ||Jj] : 
like disorder operators, they reflect the existence of non- 
trivial homotopy. From (|l]) it follows that the dimension 
of this operator is x(n; — n) = (3/2g) — 1, which is 2^2 
in the Kac classification x TyS — {{rg — s) 2 — (g — l) 2 ) /2g. 
However, for small n and n' , x(n, n') ~ (l/67r)(n — n'), so 
that, to first order, it does not matter whether x(n, — n) 
or x(0; n') is used apart from factors of 2. Note that 
if n' = —n the fusion rules of CFT imply the operator 
product expansion (OPE) a ■ a = 1 + e + — , where e 
is the (1,3) operator which has been identified as the 
local energy density of the O(n) model e(r) counts 
whether a given bond r is occupied or not - consistent 
with the interpretation that, as a — > 0, the 0(n) term in 
a(r + a)a(r — a) counts whether the single loop is trapped 
between the points r ± a. 

Here are a few examples of the application of this for- 
mula. Away from criticality it implies that (ov (i?)) ~ 
C(n; n' 

)£-x(n,n') where the 

correlation length £ behaves 
as (y c - y)- u , and (7(0; 0) = 1. When n = 0, the 
coefficient of n' counts configurations of M nested 
loops surrounding the point R of the dual lattice. De- 
noting the number of such loops of total length N by 
'■' the most singular term in the generating function 



N 

En $V ~ (1/M!)((1/8tt)| Info " v)\) M ™V^% 
using v = |. This yields the asymptotic behavior 

b^ ~ (<J /(M- l)!)(l/8^) M (l/7V)(ln7V) M -V Ar , 

where /1 = y^ 1 is the usual lattice-dependent connec- 
tive constant, and <To is a positive integer taking account 
of the fact that, on a loose-packed lattice, there are <To 
equivalent singularities on the circle \y\ = y c . For M = 1, 
one may sum over the position of the marked point R, 
instead of summing over the position of the loop, so that 
bjy = AnPn, where pn is the number of single loops per 
lattice site and is their average area. The amplitude 
1/87T in this case agrees with an earlier result found by 
a different (more involved) method |13|. It is interest- 
ing also to calculate higher point functions. For example 



the 0(n) term in the connected 4-point function (aaaa) c 
counts the number of single loops which have non-trivial 
winding around all four points. Examples are shown in 
Fig. |l|. Since a is a (1,2) operator this 4-point function 
may be computed exactly at criticality, in terms of hyper- 
geometric functions. The details will be given elsewhere 

& 

Percolation. The Coulomb gas formulation of the Q-state 
Potts model is similar to the above, except that it is valid 
only at the critical point JjJ. The model is first mapped 
onto the random cluster model, in which every cluster 
configuration is weighted by a factor of p/(l — p) for each 
bond and Q for each cluster. Each cluster may be iden- 
tified by its outer and inner hulls, which form a dense 
set of closed loops. At criticality, each hull then car- 
ries a weight \[Q. Once again this may be mapped onto 
a gas of oriented loops with phase factors e ±7rx , where 
\/Q — 2cos7rx, and then to a free field theory, where, 
however, this time g = 1 — \- The central charge then 
vanishes for g — |, corresponding to x = 3 or Q = 1; 
the percolation limit. Once again one may define a non- 
local operator which counts the hulls which surround a 
marked point with a different weight VQ 7 = 2cos7rx', 
and standard Coulomb gas methods then lead to a result 
identical in form to (|l|) for its dimension x(Q;Q'). For 
this to be given by x (1:2) = (3g/2) - 1, Q' = (2 - Q) 2 . 
For Q = 3, for example, the clusters which wrap around 
the marked point are counted with weight Q' = 1. This 
is consistent |Tj| with the identification of the Z2 twist 
operator for the 3-state model made in Ref. For 
the Ising model (Q — 2), it is the same as the disor- 
der operator. Note that as Q — » 1, the 0(Q — 1) term 
in a(r + a) o~(r — a) counts the number of clusters sep- 
arating r ± a, as compared with the 0(n) model where 
it simply counted whether or not the single loop passed 
between them. For this reason the lattice interpretation 
of the (1,3) operator in the OPE a ■ a is problematic in 
this case |la| . 

For Q = 1 + S and Q' = 1 + 8' {\S\,\5'\ < 1) 

note that x(Q;Q') - (l/4\/37r)(<5 - 5') (with 8' 28 

corresponding to xi, 2 ), while for small Q' x(l;Q') = 
Js ~ v / Q 7 + 0(Q'). As before, one application of these 
results is to count the number of clusters which surround 
a given point, that is, which have to be crossed to escape 
to infinity. If P(M) is the probability of there being ex- 
actly M such clusters, 



Y,P(M)Q' M ~C(Q')(p c -p) 



ux{l;Q') 



M 

where now v = |. For Q' = this gives -P(O), which is 
the probability that the dual site R is connected to the 
boundary by a set of dual bonds: this gives the usual 
exponent (3 = -Jr of percolation. Expanding in powers of 



ity now gives 



P(M) ~ P(0)(|ln(p c - P )\/2k) 2M /{2M)\ 
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Note that although the amplitude in P(0) is not univer- 
sal, the ratios P{M)/P(0) are. 

However, this is valid only for M <C | \n(p c — p)\. To 
find the behavior near the average value, expand around 
Q' = 1. The mean value M ~ (1/3-^/3^) | ln(p c — p)|, and 
the variance M 2 — M oc | ln(p c — p)| also. This implies 
that, as p — > p c — , the distribution of M becomes peaked 
about M. Alternatively, one could work at the percola- 
tion threshold in a finite system of size L, in which case 
M ~ (l/4\/37r) lni. 

5pm quantum Hall transition. Recently Gruzberg et al. 
pof have shown that certain properties of a model of non- 
interacting quasiparticles for the spin Hall transition (a 
2D metal-insulator transition in a disordered system in 
which time-reversal symmetry is broken but SU(2) spin- 
rotation symmetry is not pi) may be mapped exactly 
onto percolation. In particular, the mean conductance 
between two extended contacts on the boundary of a fi- 
nite system is (apart from a factor of 2 for the spin sum) 
equal to the mean number of distinct clusters whose outer 
hulls connect the two contacts. As argued above, such 
quantities are related to the derivative with respect to Q 
at Q = 1 of correlation functions of a twist operator. This 
will however now be a boundary twist operator. While it 
is possible to adapt the above Coulomb gas arguments to 
account for the boundary, since in other examples such 
methods are known to fail for boundary operators, we 
shall instead give a more direct argument. 

First consider the example of an annulus of width L 
and circumference W. The geometry is shown in Fig. ||. 
Apart from the clusters whose outer hulls cross the sam- 
ple, there are those which touch the lower edge but not 
the upper, those which do the opposite, and those which 
touch neither edge. There may also be one cluster which 
crosses the sample but which also wraps around the an- 
nulus, so that its outer hulls do not connect the contacts. 
Denote the numbers of such clusters in a given configura- 
tion of the random cluster version of the Potts model by 
N c , Ni, N 2 , N b and N w respectively. (Note that iV c = 
if N w = 1). Let Zij(Q) denote the Potts model parti- 
tion function with boundary condition of type i on the 
lower edge and j on the upper edge. The cases of inter- 
est are where i or j correspond to either free boundary 
conditions, denoted by /, or to fixed, in which the Potts 
spins on the boundary are frozen into a given state, say 
1. Then 



is the lowest scaling dimension out of all the conformal 
blocks which can propagate around the annulus with the 
given boundary conditions. When i = j this corresponds 
to the identity operator, so that Ajj = 0, but for the 
mixed case (ij) — (/l) it corresponds to the (1,2) Kac 
operator, so that A/i = Ai )2 = (Q) in the previ- 
ous notation. This identification was previously used at 
Q = 1 in Rcf. Q to compute crossing probabilities, i.e. 
the probability that N c > 0, in simply connected regions. 
Substituting into @ gives (N c ) ~ 2irA' lj2 (l)(W/ L) as 
W/L — > 00, since N w < 1. From this follows the univer- 
sal critical conductivity v / 3/2. 

At finite W/ L the corrections to the mean conductance 
are expected to be of the order of e~ n ^ 2 ' 2W / L wnere 
A 2j2 = I at Q = 1, but the full dependence requires 
knowledge of the entire operator content of the model 
for the different boundary conditions. This seems to be 
beyond the reach of current methods. However, for a sim- 
ply connected finite sample the arguments of Ref. B may 
be adapted. Consider a simply connected region with 
contacts C\C 2 and C3C4 on its boundary, as shown in 
Fig. IH The remainder of the boundary has hard wall con- 
ditions on the quasiparticle wave functions, correspond- 
ing to free boundary conditions on the Potts spins. The 
mean number of clusters crossing between the contacts 
is still given by (||) (with N w — 0), where the different 
boundary conditions are placed on the segments C\C 2 
or C3C4, with the remaining boundary Potts spins being 
free. This may then be written in terms of correlation 
functions of boundary condition changing operators |L7| 



W/i(Wi/(tf2)>(0/i(C3)&/(<7 4 )) 



These correlation functions are computed by conformally 
mapping the interior of the region to the upper half 
plane. Any conformal rescaling factors for Q 7^ 1 cancel 
in the ratio, which then depends only on the cross-ratio 
7] = [z\ — z 2 )(z 3 — Z4)/ \z\ — z?,){z2 — Z4) of the images z l of 
the points Cj under this mapping. For a rectangle with 
|CiC 2 | = W and |C 2 C 3 | = L, n = (1- k) 2 /(l + k) 2 where 
W/L = K(l - k 2 )/2K(k 2 ) and K is the complete elliptic 
integral of the first kind. Since 4>fi is a degenerate (1,2) 
operator, its four-point function satisfies a hypergeomet- 
ric equation. The details of this calculation will be given 
elsewhere uM . The result for the mean conductance is 



Zff 
Zif 



(QN C +N ul + N 1 + N 2 + N b ^ 
{Q N 2+ N b) z 



z n = (Q 

(Q N *+ N ») 




(|)m m 



so that 



(d/dQ)\ Q =i(Z ff Zn/Z n Z lf ) 



(2) 



According to the theory of boundary CFT [j^j, ~ 
exp [tt((c/24) - Aij)(W/L)] as W/L -> 00, where Ay 



For W/L 1 this reproduces the above result for the 
conductivity. In the opposite limit g ~ ^4 e ^( 7r / 3 )( i / w/ ) ; 
in agreement with Ref. J1Q| , but now with a definite pref- 
actor A = 3r(|)/2r(|)V 
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FIG. 1. Examples of loops which wind non-trivially around 4 points, counted by the connected 4-point function {aaaa) c 




FIG. 2. Annular geometry with contacts along either edge. Periodic boundary conditions are implied in the horizontal 
direction. Examples are shown of Potts clusters of types c, 1, 2 and b. The mean conductance is proportional to the mean 
number of type c. 
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c 2 

FIG. 3. Simply connected region with contacts C1C2 and C3C4 along its edge. 
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